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Recently the problem of Unambiguous State Discrimination (USD) of mixed quantum states has 
attracted much attention. So far, bounds on the optimum success probability have been derived 
1]. For two mixed states they are given in terms of the fidelity. Here we give tighter bounds as 
well as necessary and sufficient conditions for two mixed states to reach these bounds. Moreover we 
construct the corresponding optimal measurement strategies. With this result, we provide analytical 
solutions for unambiguous discrimination of a class of generic mixed states. This goes beyond known 
results which are all reducible to some pure state case. Additionally, we show that examples exist 
where the bounds cannot be reached. 



I. INTRODUCTION 



Quantum state discrimination 2] is a fundamental task in quantum information theory, especially in a communica- 
tion context. Whenever the signal states are nonorthogonal, perfect discrimination becomes impossible. One has then 
to resort to optimum state discrimination strategies by specifying figures of merit that define some optimal strategies. 
The optimum strategy depends then on the quantum states and their a priori probabilities. One strategy is Minimum 
Error Discrimination (MED) 0] in which the measurement identifies the possible input states with some error. It 
is the goal to minimize the error. Another strategy is to optimize the mutual information between the sender and 
receiver. 

The scenario studied here is Unambiguous State Discrimination (USD) which characterizes a measurement which 
either identifies a signal state without error ('unambiguous') or sends out a flag stating that it failed to identify the 
state. The objective is to minimize this failure probability. The problem of finding optimal USD strategies has been 
solved for many pure state scenarios 0, EL IE IE 1 including any two pure states [1] . 

In contrast to the MED problem, which is already solved for any pair of mixed states , optimal USD of mixed 
states is an open problem. Some special cases have be en g iven for which the corresponding problem can be reduced 
to USD of pure state case, such as in state filtering [|| flOt lllj or state comparison [J, [jjj . The underlying reduction 
theorems have been stated in [l^. F° r the general case, necessary and sufficient conditions for the optimality of a 
POVM were derived in [l3L[T^ |. They allow a numerical treatment of the problem but have not given rise to analytic 
solutions. 

For the unambiguous discrimination of a pair of mixed states, lower bounds on the failure probability have been 
found P, Il5j and reveal three regimes, depending on the ratio between the two a priori probabilities of the two 
mixed states. The boundaries of the middle regime were recently refined in |ri| but the consequences for the two 
remaining outer regimes were not addressed. Here we provide new bounds in those two regimes. Furthermore we 
derive necessary and sufficient conditions to reach the three bounds in the three different regimes. Given two density 
matrices po and p\ and their a priori probabilities rjo and 771, the necessary and sufficient conditions to reach the 
bounds given here take the form of the positivity of two particular operators. Moreover we show that examples exist 
where the bounds cannot be attained. When the necessary and sufficient conditions are fulfilled, we give the optimal 
measurement strategy to reach the bounds. 

The structure of this paper is the following. In the Sec. [HJ we derive lower bounds for the success probability in 
the case of two mixed states. Our derivation uses the Cauchy-Schwarz inequality, as used in [Tij . and allows us to 
look for necessary and sufficient conditions to reach the lower bound in each regime of the a priori probabilities. In 
Sec. IIIII we report the notion of parallel addition that leads to some useful relations for USD in connection with a 
reduction theorem of Ref.^2]- I n Sec II VI we derive the main result of this paper as a theorem: two necessary and 
sufficient conditions for the failure probability to reach the bounds are given. We also give the corresponding optimal 
POVM. In Sec0 we provide examples showing that there are generic mixed states of interest for which the necessary 
and sufficient conditions are fulfilled and for which we can therefor give the optimal USD measurement. 
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II. LOWER BOUNDS ON THE FAILURE PROBABILITY 

In Unambiguous State Discrimination, the performed measurement either identifies uniquely a state (conclusive 
result) or fails to identity it (inconclusive result). The goal is to optimize that strategy by finding the measurement 
for which the probability of inconclusive result is as small as possible. The problem is then specified by the set 
of quantum signal states {pi} and their respective a priori probabilities {r/i}. The measurement is a generalized 
measurement i.e. a Positive Operator- Valued Measure (POVM) A POVM is a set of hcrmitian and positive 
semi-definite operators {Ei}i that add up to identity acting on the Hilbert space spanned by the signal states, i.e. 
"Y^iEi = Ift. Given N possible input states, we consider measurements with N + 1 outcomes where the first N 
outcomes identify a state and the last one corresponds to inconclusive results where the identification failed. The 
POVM elements are denoted by Ek with fc = 1, . . . , N and Ei respectively. The probability to obtain an outcome for 
some POVM element E for a given signal p is then given by Tr(Ep). 

In general, a POVM describing a USD measurement satisfies Ti(EkPi) = whenever k ^ i so that only the state 
Pi can trigger the measurement outcome connected to Ek ■ The failure probability Q of a USD strategy is then given 
by Q — J2iQi> where- Qi = ^2iViTr(E?pi). From this definition we find that Qi < rji. In this paper, we consider 
the USD of two signal states po and p\ that are mixed states with a priori probabilities ijo and r/i . Accordingly, our 
POVM contains three elements {Eq, E\, E?} which correspond respectively to the conclusive detection of p , to the 
conclusive detection of p\ and to an inconclusive result. The failure probability then equals Q = Qq + Q\. 

Our interest is first focused on the product QqQ\. We can give a lower bound that is expressed in terms of the 
fidelity F of the two states. The fidelity is defined as F = ^(^^/popi^/po) The bounds, formulated in the 

following theorem, are tighter than the one given in Moreover, we pay additionally attention to the condition 
under which the bound can be reached. 

As for the notation, consider an hermitian and positive semi-definite operators O. We can define its unique square 
root \fO and decompose it into the form O = MM^ with M = VOU, for any unitary matrix U. Since the states pi 
and the POVM elements Ek all are hermitian and positive semi-definite operators, we can introduce their square root 
and use the previous decomposition. 

Theorem 1 Let p$ and p\ be two density matrices with a priori probabilities r/o and r/i . We define the fidelity of the 
two states po and pi as F = Tr{ y/^popi yfpo) ■ Then, for any USD measurement, the product of the two probabilities 
Qq and Qi to fail to identify respectively the state po and p\ is such that 

QoQi > VQViF 2 - (1) 
The equality holds if and only if the unitary operator V arising from a polar decomposition 



VPoVPi = y \fPoPiVPo V 

satisfies 

V i y/toEh=a y /n3i (2) 

for some a € M + . 

Before we turn to the proof of this theorem note that relation J2J implies a condition required for the optimality of a 
USD POVM (see [Tolll^ ). It is clear that optimality of a specific USD measurement implies that the conditional states 
after the inconclusive results do not allow further USD measurements. That would already be satisfied if, for example, 
the supports of the conditional states coincide. We find a stronger property whenever equality holds in Theorem 1. 
Indeed, if we have \fp§E-t — a^/plE-r with a € M + , then it follows immediately that \[WtpQ\f~E> = a 2 y/Ej pis/Ef , 
This means that the conditional states corresponding to inconclusive results must be identical up to normalization. 
Therefor no information whatsoever about the signal state can be extracted from these conditional states. 

Proof of Theorem 1 The basic ingredient for the derivation of the bound is the Cauchy-Schwarz inequality: 

Theorem 2 Cauchy-Schwarz inequality 

If x and y are members of a unitary space then ||a;||||?/|| > \{x, y)\. The equality holds if and only if x — ay for some 
a in C 

A unitary space is a complex linear space S together with an inner product from S x S to C. Therefore the complex 
space of bounded operators acting on a Hilbert space is a complete unitary space if we consider for two elements A, B 
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the inner product Tr{AB^). Hence, with E{ = M{M\, po — -JpEUU^ y/po and p\ = y/pi-Jpx, we obtain 

^Ti{E. lPo )^v{E lPl ) = \J Tr(U ^M t M} vW f ) \J Tr ( ^fpiMi M.\ y/pi) 
> \Tr(UJpSMiMljpd\, 

where we have used the freedom in the decomposition of po- By Theorem 2, the equality holds if and only if 

Uy/poMi = a^fpiMi, for some a G C or, equivalently, if and only if U ^fpaEi = a^/~p\Ei, for some aeC. 

We now consider a USD POVM {-Ei}i=o,i,?- Using the fact that Tr(£: pi) = Tr^po) = 0, we find for E and E x 

= ^/TiiE Q p Q )^/TiiE Q p 1 ) > \Tt{U^EoVpT)\, 
= ^Tr(E lP0 )y/Tr(EiPi) > \Tv{U^E 1 
This simply means that Ti(Uy/~poEoy/pi) = Tr(U^/poEiyfpi) = 0. For E-?, we obtain 

y/Tr(E Wo )y/Tr(E Wl ) > |Tr (U^E ? ^pT)\. 

From this it follows that we can write 

^Tr(E ?Po WTr(E ?Pl ) > \Tv(U^E,^) + + 0| = \Tr(U^^)\ , (3) 

where we used the relation yV Ej = 1. Furthermore, the inequality Q must hold for any unitary matrix U so that 
we find 

^Tr(E Wo )^/TiiE ? p 1 ) > max|Tr(C/V^VPT)l- ( 4 ) 
Here, again, the equality holds if and only if a unitary operator J7 max which maximizes the right hand side satisfies 

tAnax\/P0-E? = OL^pxEj 

for some a G C. To find the unitary matrices J7 max that maximize \Tr(U, l /poy/pi)\ we use the following lemma: 
Lemma 1 For any operator A in the space M n of n x n matrices we find 

m&x\Tr(AU)\ = Tr(\A\) 

where the maximum is taken over all the unitary matrices. The maximum is reached for any unitary operator U that 
can be written as U — V^e 1 ^. Here is an arbitrary phase while the unitary operator V is defined via the polar 
decomposition 

A=\A\V 

with \A\ = \f~AA^ = vVAfAV^. 

Proof For any operator A, we can introduce its polar decomposition A = \A\ V with \ A\ — VAA^ = V VaTA . 
Note that V is unitary while V AA^ and V A^A are unique, positive semi-definite and hermitian. With that we find 

\Tr(AU)\ = |Tr(|j4|FE/)| = |Tr(|v4| 1/2 |^l| 1/2 U[/)|. 

We denote X = \A\ 1 / 2 = X^ and Y = \A\ 1 / 2 VU and apply the Cauchy-Schwarz inequality (Theorem 2) to obtain 

\Tr(AU)\ = |TY(Xty)| < V^W) \]Ty{UW\A\VU)) = Tr(|A|) . 

Equality holds if and only if IA] 1 / 2 = /3\A\^ 2 VU, for some /3 G C. This is possible if and only if (3VU = 1, where U 
and V are both unitary matrices. This means that f3 = er 11 ^ for some <fi so that we find the connection U — V^e 1 ^ . 
This completes the proof of the lemma. 
Thanks to lemma ^ Eqn. 0} implies 

v/Tr(i? ? po)VTr(i? ? pi) > ^(IV^V^D! 
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where equality now holds if and only if 

Vie*y/nEh = a^Ei (5) 

for some a £ C. Let us introduce the operators Fq := \^fpo^/pi\ = \J \[PoP\ \[Po an d F% = V^FqV — J \fp~iPoyfpi, 
which are motivated by the polar decomposition 

pE^pl = F Q V = VFx. (6) 



These operators are related to the fidelity of the two density matrices through the relation F = Tr(y^/pTpo^7pT)(= 
Tr(Fo) = Tr(Fi)) ([HI). 

Next we use the definitions of the partial failure probabilities Qi — ?/;Tr (E->pi) and choose the phase e 1 ^ to be the 
same as the phase of a in JSJ to obtain the desired inequality QoQi > rjorjiF 2 . Equality in the previous equation 
then holds if and only if V^i/poE-? — u^fp\E>, for some a e M + . This completes the proof. ■ 

We can now derive the bounds in the different regimes of the ratio between the two a priori probabilities. 
Actually, the procedure is to find the minimum of the failure probability Q = Qo + Qi under the constraints of the 
previous derived inequality QoQi > rjorjiF 2 . According to Theorem 1, we can provide the necessary and sufficient 
condition for equality. 

Theorem 3 Let p and pi be two density matrices with a priori probabilities r/o and r\\. We define the fidelity F of 
the two states po and p\ as Tr( *J ^/popi yfpo) ■ We denote by Pq and Pi, the projectors onto the support of po and p%. 
Then, for any USD measurement, the failure probability Q obeys 

Tr(P lPo ) V rj F 

Q>2y/Wh F for & <\—< 



F ~ V Vo ~ Tr(Popi) 



F 2 F Pr\f 

Q > m^rrB \ + Vi Tr{P a pi) for — — < J — 

Tr(Popi) Tr[Po P i) ]/ r) 



Equality holds if and only if the unitary operator V arising from a polar decomposition ^/pa^/pi = \J \fpo~P\\/~Po V 
satisfies \fp§Ei = a^/piE-?, with a = ]jlE±£lL ^ a — FEL an d a — F — - { n the the first, second and third regime, 



lo Tr(Papi) 

respectively. 

Proof First of all, according to Theorem 1, we know that for any USD measurement the inequality Q\ > ^''q^ 
holds. It follows that the failure probability is such that 

Q>Qo + ^. (8) 

Qo 

Let us consider relations that only hold if equality holds in Eqn. (JHJ). In this case we have 

QoQi = VoViF 2 • (9) 

Moreover, from Theorem 1 we know that in this case we have W \fp~Q~Et — ot^fp\E>, for some a € K + . This relationship 
implies, via the respective definitions, that 

Qo = « 2 -Qi- (10) 
m 

We can combine the two equations and (|1C)|) to 

Qo = ar/oF . (11) 

So the final statement is that Q = Qq + 7,0 q 1 ^ if and only if ^fp^E-t — a^fplE-i , where a now is explicitly related 
to the other parameters as Qo = arjoF. 

Second, we have to derive the range constraint on Q and Q\. We know already that Qi < r^. Moreover, from the 
work by Herzog and Bergou in |ll|. we learn that t]qTi(PiPo) < Qo an d ?7iTr(Popi) < Q\. Indeed, from the structure 
of the POVM elements, we have Eq + E\ + E-> = 1 with Se C and Sbx C S^ a . We consider only the non-trivial 
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case where the supports of po and p\ are not identical. Then its structure must be such that E± +E-r — P\ + R where Pi 
is the projection onto the support of pi and R is an hermitian positive semi-definite operator with support Sr C 

which satisfies E + R = P^. Then it follows that P = r] Ti(E p a ) = rj Tr(P^p ) - r) Tr(Rp ). In our non-trivial 
case we will have Tr(Ppo) > as soon as R ^ 0. This yields Po < ^oTi^P^po) or equivalently Qo > ^oTr(Pip )- In 
the same way, on can find Q\ > ?7iTr(Popi). We then have 

77oTr(P lP o) < Qo < Vo, (12) 
?7iTr(Popi) < Qi < rji. 

These two constraints can be combined to 77 Tr(Pi/9o) < Qo < Vo t f n — \ • This can be seen as follows. Since 

Qi = n ° r !i F , the constraints on Qi take the form rj F 2 < Q < r] rp-f- — -. Let us consider the USD POVM given 
" Wo " ir(Popi) 

by {E-? — P\,Eq = Py-,Ei = 0}. Thank to Theorem 1, we find r] r]iF 2 < ry ?7iTr(P 1 /9o)Tr(P 1 p 1 ) or in other words 
VoF 2 < ?7oTr(Pipo)- We can also consider the USD POVM given by {E ? = P ,E =0,E l = P^} and with Theorem 
1, we finally have ?7otv75 — T ^ Vo- 



Ir(PoPi) 

Next, we define the function q(Qo) = Qo + '"'q 1 ^ and minimize it under the constraint rjoTr(Pipo) < Qq < 



r/o rp f — -. The resulting minimum will constitute a lower bound for Q. The function q(Qo) is convex (4r&(Qo) > 

0) and, therefore, it takes its minimum at the point Q™ 111 where the derivative vanishes (^^-(Qo) = yielding 

Qo 1 ™ = VVomF) or at the limits of the constraint interval (QS lin = ?7 Tr(Pip ) and QS lin = ^ot^ — That 

ir(Popi) 

gives us the minimum in three different regimes. In the first regime we have q m in(Qo) — ^oTrfPipo) + Vi t f n — r 
and Q nm = tjqTt(Pi/9o) if y/rjo~rjiF < t]qTt(PiPq) that is to say if Jl^ < Tr(Pip n ) ^ j n ^ c second regime we have 



<7min(Qo) = 2^jo~mF and Q$ in = ^mmF if Tr( ^ lpo) < < Tr( p opi) - The third regime gives 9min (Q ) = 

"Tffc + ^{Jh/n) and (tf* = Vot^ if < Jf - 

As a result we obtain lower bounds for the failure probability Q in three regimes as given in Eqn. J7J. For each 
regime, the value of Qo which minimized q(Qo) is given and via Eqn. i|llf) we find the corresponding value that a has 

to take. We read off the values as a = Tr(Pipo) _ /?jt anc j Q _ F — f or the nrs t second and third regime, 

f ' y io Ir(Popi) 
respectively. ■ 
Let us note that, by construction, those bounds are tighter than the ones in Q]. Indeed, one could recover the three 
bounds in by looking for the minimum of the function q(Qo) under the weaker constraints 770P 2 < Qo < Vo- 



III. THE PARALLEL ADDITION p Q : p x 

Before deriving our central theorem, we will first recall some useful results of linear algebra. We denote by M _1 
the pseudo-inverse of a matrix M, which has not necessarily full rank. The pseudo-inverse can be defined via the 
singular-value decomposition of M. Whenever M is of full rank, the pseudo-inverse coincides with the inverse. In 
general, it is not known how to express the pseudo inverse of a sum (A + B)^ 1 in terms of the pseudo inverses A^ 1 
and B^ 1 19]. However, a related new operation A(A + B)~ 1 B 1 called parallel addition and denoted by A : B has 
been defined in 1969 by Anderson and Duffin and will turn out useful in our context. 

First of all, we denote by Sm, the support of a hermitian and positive semi-definite matrix M. We then have the 
following property for the parallel addition: 

Property 1 poj Let A and B be two hermitian and positive semi-definite matrices in M n , then the support Sa-.b of 
A: B is given in terms of the supports of A and B as 

Sa-.b = Sa^Sb- 

Next let us recall two reduction theorem for USD of mixed states 0] . We consider the problem of discriminating 
unambiguously two density matrices po and p\ with a priori probabilities rjo and rji . We denote by rg the rank of po 
and by r\ the rank of p\ . A general USD problem can satisfy ro + t\ > d, where d is the dimension of the Hilbert 
space 7i spanned by the two states. This means in particular that the two supports can overlap. 

In a first reduction theorem it has been shown by the authors 0] that any such USD problem can always be 
reduced to the one of discriminating p' and p[ , two density matrices of rank r and r[ with a priori probabilities ry 
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and rf-y , spanning the same Hilbert space TL of dimension r' Q + r[ . Indeed we can split off any common subspace of 
the supports S po H S pi to end up with S p > o H 5^/ = {0}. An easy way to know whether the two supports overlap is to 
check whether the equality rk(p' ) +rk(p' 1 ) = rk(p' + p[) holds (see details in [1]]]). In the reduced case, property 
implies Sp'g-.p^ — that is to say p' : p\ = 0. By defining £ := p' + p' x , we can write the parallel addition as p' Y,~ 1 p 1 . 
Moreover, since rk(p' + p[) = dim(7i), we end up with £ having full-rank and ££ _1 = 1^. 
We therefore have the following corollary to property (1), 

Corollary 1 Let po and p\ be two density matrices spanning a Hilbert space Ti.. Let £ be defined as the sum of these 
two density matrices. 

If rk(po) + rk(px) = rk(p + pi) then po^ 1 Pi = 0. 

According to the first reduction theorem we can, without loss of generality, consider only USD problems of two 
density matrices without overlap of their supports. In the following, we consider two density matrices po and p\ 
(which are hermitian and positive semi-definite matrices) such that rk{po + p{) = rk(po) + rk(pi) = dim(7i). As 
explained above, for such a problem, po^Vi = 0, with £ = po + p\ having full rank. This leads to po'E~ 1 po = Pa 
and /9i£ _1 pi = pi since ££ _1 = 1^. The projectors onto the supports of those two density matrices can then be 
written as : P Pa = ,/ft)£~\/pn and P Pl = ^jTT 1 ^fp~[. 

A second reduction theorem [l^] allows to eliminate the part of the support of p' which is orthogonal to the support 
of p[ and vice et versa. This implies that the two resulting density matrices p' ' and p'[ possess the same rank r and 
span a 2r-dimensional Hilbert space. We denote such a USD problem by "r + r = 2r" (see for more details). This 
second reduction theorem indicates in which situations a further reduction of the original problem can be achieved. 
This theorem is not needed for the derivation of our central theorem. 



IV. NECESSARY AND SUFFICIENT CONDITIONS 



We are now ready to derive the main result of this paper. The first part of this result gives compact necessary and 
sufficient conditions for a pair of mixed states to saturate the bounds of the failure probability Q. The second part 
gives the corresponding POVMs in an explicit form. To clarify the notation, let us note that in M > we mean that 
the operator M is hermitian and positive semi-definite. 

Theorem 4 Necessary and sufficient conditions to saturate the bounds on the failure probability 
Consider a USD problem defined by the two density matrices p Q and p\ and their respective a priori probabilities r)o 
and r]i such that their supports satisfy S Po D S Pl = {0} (Any USD problem of two density matrices can be reduced to 
such a form according to fZl/). Let Fo and F\ be the two operators y/ \[PqP\\[Pq and y/y/piPo^Pi- The fidelity F of 
the two states po and p\ is then given by Tr(Fo) = Tr(F\). We denote by Pq and Pi, the projectors onto the support 
of po and p\. The optimal failure probability Q°P^ for USD then satisfies 



Q 



opt 



Vi T( F p +mTr(P lPo ) P n ° 
Tr(P 1 p ) pi 



Tr(P lP0 



Tr( Pl 



, ^o>0 , fm < Tr(P lPa ) 
* *i>0 f ° r V^-^^ 



Q opt = 2^/rjol^F «■ 



Po 



Pi 



^-F > 

no 

m 1 



for 



Tr(P lPo ) 
F 



< 



< 



F 



T]o Tr(P pi) 



(13) 
(14) 



Q 



opt 



Tr(P pi 



b 

+ m Tr(P QPl ) * P °~ ggjjg 

F 



Pi 



1 Fo>0 
Fi>0 



for 



F 



Tr(PoPi) 



< 



The POVM elements that realize these optimal failure probabilities, if the corresponding conditions are fulfilled, are 
given by 



E = £ 1 Vpo (fio - aFo) % /p^£ 1 



E 1 = S^VpI ( Pi ~ -Fx ) VpIS" 



(15) 



with a = FfiPiPa) j or ^ e ji rs i re gi m e, a 



for the second regime and a — jy^p — r f or the third regime. 
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Proof First, we give a proof for the necessary conditions. 

Proof for the necessary conditions From Theorem |21 we know that the bounds on the failure probability are 
satisfied whenever VK/poE-? — a^/pTE? with a = TEi^fsl a = IHL anc [ a — F — f or the three regimes. 
We replace E? by 1 — Eq — Ex, multiply on the left by V and on the right by Jp~^. This leads us to 



p - aF = VPo-^oa/Po (16) 

where we used the relation (J5J) ^/po^/pi — FqV and the fact that the support of pi and Ej are orthogonal for i ^ j. 
Indeed, let us notice that Tv(EiPj) = EiPj = because B, and pj are hermitian and positive semi-definite 
operators The right hand side in 1)16(1 is hermitian and positive semi-definite because of the form AA^ with 

A = y^PoV -Bo- Then po — c*-Fo must be hermitian and positive semi-definite as well. A similar calculation where we 
only multiply on the right by v / pT instead of by v / po leads us to 



pi - -F x = ^fplEx^fpl 
a 

which is again a hermitian and positive semi-definite operator. 

With that we have proved that if equality holds in the bounds of Theorem [3] then we have 

Po-aF > (17) 

Pi --Fx > 0, 
a 

which form, therefore, necessary conditions for equality in the bounds of Theorem [3J 

Proof for the sufficient conditions Now we start with the assumption that the conditions 117(1 are fulfilled. Let 
us define the following POVM elements : 



E = ^-^(po-aFo)^^- 1 (18) 
Ex = Z- 1 ^px~(px--Fx) VPTST 1 



a 



First, let us verify that this is indeed a valid POVM. The three operators are positive since they are of the form 
A] MA where M is a positive hermitian operator. In the first two cases this is true because of the conditions 1(171) . in 
the third case it follows from the positivity of F . The three operators sum up to identity, E + Ex + E-> = 1, as can 
be checked by straight forward calculation which makes use also of Eqn. JfjJ. Next, we have to check that the given 
POVM is a valid USD POVM, that is, Tr(p Q Ex) = Tx(pxE ) = 0. This relation holds since the supports of p and pi 
do not overlap. Therefore, corollary 1 applies and we have poE _1 pi = from which follows that y/po^ ~Vi = an d 
y/piE _1 po = 0. Finally, one can check in a straight forward calculation exploiting the properties used in the previous 
checks that this POVM lead to the three desired failure probabilities. 

Let us note that we have only used the assumption about the non-overlapping supports to prove the sufficiency of 
the conditions. Their necessity does not require this assumption. ■ 



V. DISCUSSION 



Theorem 0] characterizes under which circumstances the equality of the bounds in Theorem [3] can be obtained. 
Whenever two mixed density matrices have no overlapping supports and the corresponding two operators in Theorem 
0]are positive semidefinite, we can give explicitly the optimum USD POVM. 

The first question is to know whether the set of pairs of generic mixed states (a USD problem which is not reducible 

to some pure state case), that fulfill the constraints po — ^J^Fq > and px — ^J~^Fx > 0, is empty or not. Actually 

this set is non-empty. For instance, consider a problem motivated by a four-state quantum key distribution protocol 
using coherent states p3 |. Here it might be of interest for an eavesdropper to distinguish the density matrices 
po = 5 [|a)(a| + | — a|] and pi — ^[\ta) (ia \ + \ — ia){— ia\], corresponding to the bit value and 1, respectively. In 
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fact, this pair of states can be represented as geometrically uniform (GU) states |2^| as they are related as p\ = UpoW 
with U 2 — 1. They can be represented as operators over a four-dimensional Hilbert space as 



Po 



( M 2 
o 



C0C2* 

ci | 2 C1C3* 

c 2 c * |c 2 | 2 

V c 3Cl * |c 3 | 2 



(19) 



with complex coefficients Ci depending on phase and as given in |22| . and 



U 



-1 \ 
0-100 

10 

1 / 



(20) 



One can show that for these two states the operators po — yT^Fo and p\ — y jj^-^i are hermitian and positive semi- 
definite for some regime of the ratio — around the value — 
C3 in C 

F = e 



1 included into the second regime (for any cq, c\, c 2 and 
According to Theorem the optimal failure probability is Q opt = 2^f\^r\{F where the fidelity is given by 



" (| cos 



+ sin 



Let us note that, in general, those operators are not positive for the whole second 



regime F < y^ < -p- Actually this depends on the parameters Cq, C\, c 2 and c 3 in C. This implies that, in general, 

the necessary and sufficient conditions are not fulfilled neither for the first regime nor for the third regime for these 
two coherent states. 

Actually two GU states are not necessarily in the set of states that saturate the bound, not even for equal a priori 
probabilities. For example, one can consider the two GU states po and p[ = WpoW^ where po is given as above while 
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(21) 



with Co = %/cTI, ci = V0.4, c 2 = V013 and c 3 = V0.2 and rjo = rj\. Those states are indeed GU states since 
W 2 — 1. The supports do not overlap. However, one can show that the operators pa — Fq and p\ — F\ are not positive 
scmi-dcfinitc. 

As a result, there exist generic mixed states that satisfy the conditions of Theorem 4 and for which a optimal USD 
strategy can be given. However, there are generic mixed states that do not satisfy the conditions so that it remains 
to find the optimal failure probability in those cases. 



The second remark is about the link between our result and the pure state case. Actually for two pure states, since 
F = F|* )(*o|, ^i = F|*i>(*i| and Tr(P Pi) = Tr(P 1 p ) = F 2 , the constraints p a - aF > 0, p x - ~F X > are 
always fulfilled and our result reduces to the one of Shimony and Jaeger. We can go beyond this remark and find under 
which conditions our bounds reduce to the ones in 0. Since our bounds are tighter, the bounds in are reached if 
and only if, first, the condition in Theorem 4 are fulfilled and, second, the equalities Tr(P Pi) = Tr(P ly o ) = F 2 hold 
(like in the pure state case). This is made more precise in the following corollary to Theorem 4: 

Corollary 2 Necessary and sufficient conditions to saturate the bounds in 0/ 

Consider a USD problem defined by the two density matrices po and p\ and their respective a priori probabilities t]q 
and rjx such that their supports satisfy S Po PI S Pl = {0} (Any USD problem of two density matrices can be reduced to 
such a form according to fZl/). Let Fq and F\ be the two operators */ ^PoPi^Po and ^y/piPo^Pi- The fidelity F of 
the two states po and p\ is then given by Tt(Fq) = Tr[F\). We denote by Pq and Pi, the projectors onto the support 
of po and pi. The optimal failure probability Q°P^ for USD then satisfies 
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opt 
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(22) 
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Pi - FFx > 



The POVM elements that realize these optimal failure probabilities, if the corresponding conditions are fulfilled, are 
given by 

E = Z-^ipo-aFo)^^ 1 (24) 
with a = F for the first regime, a = x / — for the second regime and a = i for the third regime. 



In the first regime we find that E\ = because this operator is hermitian, positive semi-definite and its trace 
vanishes. The resulting POVM has to be a projective measurement with projections onto the support of p\ and onto 
its orthogonal complement, i.e. Eq = Pf~, E% = and E-> = P%. A direct proof from the explicit expressions in Eqn. 
(11811 is difficult, however a simple reasoning allows to verify this statement. We consider only the non-trivial case where 
the supports of pi and pi are not identical. Of course, a two-element USD POVM satisfies Eq+E? = 1 with Se C S Pl . 
Then its structure must be such that E-? = P% + R where Pi is the projection onto the support of p\ and R is an 
operator with support Sr C which satisfies Eq + R = P±, Then it follows that Q = r\\ + rjoTr(Pipo) + r/oTr(Rpo). 
In our non-trivial case we will have Tr(Ppo) > as soon as R ^ 0. Therefore we find as an optimal solution within this 
class of two-element USD POVM, the POVM with R = leading to P? = Pi and E = P^. We can actually write the 
failure probability as Q opt = r)i + i]oF 2 . Indeed p\ = yF\ then p\ = ^^fpiPQ^fpi- This implies F 2 pi = PipoPi and 
finally Tr(Pipo) = F 2 . This is consistent with the results derived above and gives the correct failure probability. In 
the third regime, we have Eq = and the corresponding POVM is a projective measurement with Eq —Q,Ei= Pq-, 
E, = P Q . 

Finally, let us note that the optimal error probability for the minimum error discrimination strategy is Qmed = 
| (1 — Tr(|77ipi — r/ po\) Q{. Then, on one hand the trace distance is related to the minimum error discrimination 
while on the other hand the Fidelity is related to the unambiguous state discrimination strategy. 



VI. CONCLUSION 



To summarize, we have given new bounds on the failure probability of unambiguously discriminating two mixed 
states. Moreover, we provide necessary and sufficient conditions for two mixed states to saturate those bounds. With 
that result, we give the optimal USD POVM of a wide class of pairs of mixed states. This class corresponds to pairs 
of mixed states for which the lower bounds (one for each of the three regimes depending on the ratio between the a 
priori probabilities) on the failure probability Q are saturated. This class in non empty since it contains some pairs 
of generic mixed states as well as any pair of pure states. For those pairs, we provide the first analytical solutions for 
unambiguous discrimination of generic mixed states. This goes beyond known results which are all reducible to some 
pure state case. Additionally, we showed that there exists pairs of mixed states that cannot saturate the bounds. 
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